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We propose and analyze two series of clustered quantum Hall states for rotating systems of spin-1 
bosons. The first series (labelled SU(4)fc) includes the exact ground states of a model Hamiltonian 
at large angular momentum L, and also for N — 3k particles at L = A''. The latter is a spin-singlet 
boson triplet condensate (BTC). The second series, labelled SO(5)fc, includes exact ground states 
at large L for diff'erent parameter values. 



In the study of the phenomenon of Bose-Emstein con- 
densation (BEC) of cold atoms, the possibility of using 
bosons with internal degrees of freedom (spin 1 in par- 
ticular), and of rotating the condensate has led to fas- 
cinating new possibilities for the states of matter that 
can be formed and observed. A multicomponent BEC 
can be realized by trapping higher spin atoms such as 
®^Rb and ^^Na ^ in optical traps, which do not affect 
the spin degeneracy, and spin-1 or "vector" EEC's have 
been realized and studied In a model description, 

the 2-body interactions of the spin-1 atoms are written 
as (5(ri - r2)[co + C2S1 • S2], with co = {go + 2g2)/3, 
C2 = (52 — 5o)/3, and gs — 47r ?i^as/M{, with Mf, the bo- 
son mass and as the s-wave scattering length in the total 
spin S channel |^,||. In the case of ^'^Na one has C2 > 0, 
and hence "anti-ferromagnetic" spin correlations. The 
"polar" ground state |^,|^ has (S) = 0; it supports many 
interesting collective excitations (see ^ for a review). 
The case of *^Rb has C2 < 0, leading to a spin-polarized 
(ferromagnetic) ground state. 

The effect of rotating a "scalar" BEC of spinless or 
spin-polarized atoms in an isotropic harmonic trap of fre- 
quency luq can be studied at mean field level (the Gross- 
Pitaevskii equation), and a vortex lattice is found at suf- 
ficiently high rotation frequency uj, in agreement with ex- 
perimental observations 0. At still higher frequencies, 
quantum fiuctuations become important and it has been 
argued that the vortex lattice is replaced by a sequence of 
distinct quantum fiuids, which can be understood using 
a mapping to the fractional quantum Hall (qH) effect . 
In this paper, we explore the corresponding possibilities 
in a rotating system of spin-1 bosons. 

In the regime luq — uj <^ luq and very weak interac- 
tion, only the lowest-energy single-particle states in the 
trap are of interest. These have non- negative angular 
momentum to = 0, 1, ... about the rotation axis, and no 
excitation of the motion along the axis, so the motion is 
effectively in the xy-plane. The x, y dependence of these 
wave functions is ipmiz) oc z^e"'^' (where z = x -\-iy, 
and the quantum length (fi./Mf,wo)^^^ of the trap, and %, 
have been set to 1); they correspond to the lowest Landau 



level (LLL) g. 

We consider the interaction Hamiltonian for bosons i, 
i — 1, . . . , N, restricted to the LLL, 



Hint = ^ S{ri - rj)[co + 028^ • Sj] 



(1) 



The ground state (GS) for the system rotating at fre- 
quency Lu is found by minimizing i/tot = ^^int + i^o — ^)L 
1^, where L — total angular momentum. 

For N particles in the LLL with zero center-of-mass an- 
gular momentum, the average "filling factor" of the oc- 
cupied states can be defined as = N/{Nv + 1), where 
Ny = 2L/N is the number of vortices there would be 
if there were a Bose condensate (we assume here a uni- 
form average occupation of states up to to = Ny). The 
use of only the LLL states is physically reasonable when 
luq ~ ^ ^o/Nv (so that particles in the LLL with 
TO = Nv cannot lower their energy by moving to a non- 
LLL state with m = 0) and vcs ^ t^^o, 5 = 0, 2 (so that 
perturbative corrections from mixing non-LLL states into 
the GS are negfigible). 

For the scalar case (in the LLL, without C2), a study |^ 
using periodic boundary conditions found that the vortex 
lattice melts at a critical v — Vc ^ 10, that incompress- 
ible fiuid states occur ati/ = i,l,|,2,...<j/c, and that 
the corresponding GS wavefunctions have large overlaps 
with the Read-Rezayi (RR) |Q series of fractional qH 
states. 

Here, for the case of vector bosons, we first propose 
two series of spin-singlet generalizations of the RR states 
(labelled SU(4)fe and S0(5)fc, respectively, and = 1, 2, 
...), and explain their structure. In a thermodynamic 
limit N 00 with v fixed, they each represent incom- 
pressible liquid phases of the system, with v = 3fc/4 for 
the SU(4)fc states and v — k ior the S0(5)fc series. Then 
we present results for the GS's of Hint as a function of L, 
Co, and C2. We identify a region in the C0-C2 plane where 
the proposed qH states with k = 1 give the exact GS's 
at high values of L. We demonstrate that the k ~ N/3 
member of the SU(4)fc series is an exact eigenstate and, 
apparently, the GS of Hint with C2 = 0, cq > 0, at angular 



1 



momentum L = N. In second quantization, this state, 
a spin-singlet boson-triplet condensate (BTC), takes the 
form 



|BTC) 



(2) 



where bl^^ is a creation operator for the single-particle 
state ipmiz) with spin /i, where ii — x, y, z. Finally, we 
show the GS's of i/tot as a function of lu. 

We now explain our analysis, starting with the con- 
struction of spin-singlet qH states for spin-1 bosons. This 
construction is a generalization of similar constructions 
for scalar and spin- 1/2 particles (either of which could be 
electrons or bosons) [|o|-^ . These papers all employed 
a correspondence between qH states and specific con- 
formal field theories (CFT's). This connection allows one 
to obtain trial qH wavefunctions as chiral correlators in 
CFT's that are associated to a Lie algebra G and inte- 
gers A: > 1 and M > 0, where M is even (odd) for bosons 
(fermions). These states generalize the pairing familiar 
from the theory of superconductivity to the formation of 
"clusters" . The structure of these CFT's is such that the 
resulting qH wavefunctions have a property that guaran- 
tees that they are exact zero-energy eigenstates of certain 
model Hamiltonians. For example, starting from G = 
SU(n -f 1) and putting M = 0, one can find completely 
symmetric wavefunctions that obey 



■^{zi,Z2, . . . , zn) ^ for zi = 
^'(zi, 22, . . . , zat) = for = 



= zi{l< fc), 

= Zk+i , (3) 



independent of the spins of the particles involved. Hence, 
these are zero-energy eigenstates of a Hamiltonian with a 
repulsive k + 1-body i5-function interaction pO| ] , and are 
the unique states of the lowest L with this property. The 
existence of such a Hamiltonian and the fact that it has a 
gap in its energy spectrum ensures that a corresponding 
incompressible liquid phase of matter exists (when N — ^ 
oo at fixed k, M) over a range of interaction parameters, 
not just for the k + 1-body model. The filling factors of 
these SU(n + l)fc states are 



u{n,k,M) 



ik 



nkM + n+l 



(4) 



The parameter n corresponds to the number of compo- 
nents of an internal degree of freedom of the particles, 
and these models have U(n) symmetry; the clusters con- 
tain nk particles. For n = I this construction gives the 
RR states, while for n = 2 it produces spin-singiet states 
for spin- 1/2 particles ||l^. Finally, clustered qH states 
admit excitations (quasiparticles) of effectively fractional 
"charge" (i.e. particle number). The simplest of these can 
be viewed as the result of the adiabatic insertion of a frac- 
tion of a quantum of magnetic flux (or vorticity), where 
the allowed fraction is 1/nk in the present cases. States 
with more than three quasiparticles at well-separated 



fixed positions display large degeneracies, which give rise 
to "non-abelian statistics" ||l^ and are understood for 



1, 2 111 



Our first series of clustered spin-singlet qH states of 
spin-1 bosons consists of the SU(4)/j states (n = 3 above), 
where we put M = 0. The filling factor is v — 3fc/4. The 
CFT construction guarantees that, for TV divisible by 3k, 
this state is a singlet under SU(3), and hence also under 
its S0(3) subgroup, the usual spin rotation group. The 
wavefunction for N — 3fcp particles can be written in 
terms of its components for particular spin states fj, = x, 
y, z specified for each particle (similar to ||lHl), as 



SU(4) 



groups 



Pgr 



groups 



2,2,2,1.1,1 



(5) 



where 



$2,2,2,1,1,1^^.^ 

n 
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^Mn2 



n 



(6) 



In words, the operations Pgroups and iSgroups in eq. (|^) 
tell one to divide the 3pfc bosons into k groups of ip 
bosons each {p of each spin polarization), to write a fac- 
tor v[/2,2,2,i,i,i gj^f-jj group, and finally, 
to symmetrize over all ways the particles can be divided 
over the k groups. For fc = 1 there is a single group, 
and one finds a state of total degree L = 3p(2p — 1) that 
generalizes the Laughlin and Halperin states for scalar 
and spin-1/2 particles, respectively, and which is a zero- 
energy eigenstatc of i?int for all cq, C2, and the unique GS 
when both ijo, 52 > since iJint is then positive. Putting 
instead k = N/3 gives k groups of 3 particles each, the 
resulting state being the P = iV BTC, eq. @. The fun- 
damental quasiparticles over the SU(4)fc qH liquids have 
fractional "charge" equal to ±1/4, and spin 1, for all k. 
We note however that such assignments, while meaning- 
ful for TV ^ fc, may be meaningless when N ^ 3k, where 
the quasiparticle size is comparable to the size of the fluid 
drop. 

An important point is that for C2 — 0, the Hamiltonian 
(|l|) has SU(3), not just SO (3), spin- rotation symmetry. 
Therefore, we expect the SU(4)fc states to be relevant 
near this line when cq > 0, for sufficiently large L. 

Our second series of qH states for spin-1 bosons is ob- 
tained from a CFT with G = SO (5), and generalizes 
the construction in Ref. [|2|. The GS's possess a sym- 
metry under an S0(3) subgroup of S0(5). The general 
construction, which involves spin-singlet clusters of 2fc 
particles, gives states of filling factor v = k/{kM + 1). 
Putting M = 0, we obtain qH states with v — k. For 
fc = 1, the wavefunction with values in the spin space of 
N spin-1 particles can be written as 



^lT{z.,)^Vi 



PiPj + <7i(Tj 



■ TiTi 



(Z, - Zj) 



■ (7) 
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Here ^'^(zi) = Y[i<ji'^i~'^j) ^ spin-independent Laugh- 
lin factor in all N particle coordinates; the Pfaffian Pf is 
defined, as usual, by 



Pf (My) = ^(Mi2Af34...M. 



N-l.N 



(8) 



where Mij are the elements of an antisymmetric matrix, 
and A denotes the operation of antisymmetrization; p^, 
(Ti, Ti are basis vectors in the spin space for particle i, 
that correspond to x, y, z, and the product is the tensor 
product. For k — 1, N must be even, and the cluster- 
ing (or pairing) of particles is seen explicitly. Clearly, 
there must be even (but otherwise arbitrary) numbers of 
particles in each of the three spin states p, a, t. The 
k = 1 wavefunction for N ~ 2q particles has total degree 
L = 2q{q - 1). The function (0) can be viewed as spin- 
singlet p-wave pairing of composite fermions of spin 1 
(for a review, see Ref. ||l^). There are also excited states 
with unpaired "charge" -neutral spin-1 fermions. We note 
that, in the state (Q), two bosons are found at the same 
point only if the total spin of the pair is zero, not if it is 
2, and hence the state is an exact zero-energy eigenstate, 
of lowest L, for a (S-function interaction that includes a 
projection onto spin 2. In our parameterization, that is 
C2 = Co/2 or go — 0, and if also Cq > 0, then Hmt is posi- 
tive, so our k = 1 state is the GS at this L. For general 
k, the S0(5)fc wavefunction can be written as a CFT chi- 
ral correlator, and is an exact zero-energy eigenstate of 
a certain k + 1-body (5-function interaction (details will 
appear elsewhere). In general, the quasiparticles over the 
S0(5)fe qH state have "charge" ±1/2 and spin 1/2, thus 
displaying a fractionalization of both "charge" and spin. 

We next report on our study (numerical and analytic) 
of the GS phase diagram of the model Hmt- The numer- 
ical work is restricted to small values of N (up to 12), 
but it indicates many features which we believe to hold 
for general values of TV. To guide our discussion we have 
displayed in Fig. |l| various special directions (rays) and 
regions in the C0-C2 plane. 

For L = 0, the GS has total spin S" = TV for C2 < 
(ferro regime) and S = (1) for N even (resp., odd) for 
C2 > (anti- ferro regime). For C2 = 0, there is a single 
SU(3) multiplet of spin states, decomposing into unique 
SO (3) multiplets of each spin S = N, N ~2, .... 

As L increases, these two phases survive in part of the 
phase diagram, as compact drops of fluid, with the center 
of mass carrying all the angular momentum. Meanwhile, 
the positive cq axis gradually opens into a region that 
contains other phases. By the time L is > A^, the co~C2 
plane contains the three regions labelled la, lb, and II 
m Fig. |. The GS's in regions la and lb are similar to 
the GS in the "attractive" regime in the scalar case 
The orbital part of the GS wavefunction is of the form 
^'(zi) cx z^, with Zc = Zi/N the center of mass (CM) 
coordinate. In region la (cq < 0, C2 > 0), the spin state 
is the same spin-singlet as for the L = GS, and the 
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FIG. 1. Overview of C0-C2 plane, with special regions and 
directions marked. 



GS energy becomes {cqN{N - 1)/2 - iVc2](27r)-3/2 0. 
In region lb (ci < 0, co < — C2), the spin state is ferro, 
S" = TV, giving energy (co + C2)N{N - l)(27r)-3/2/2 < 
At C2 = 0, Co < 0, the spin states again form the 
SU(3) multiplet. In the "repulsive" region II, the GS 
is in general not a common eigenstate of the co and C2 
parts of the interaction, and the GS energy depends non- 
linearly on the ratio cijcQ. The L — N GS's all have 
5 = or 5 = 1 (depending on the value of A^ modulo 
6, and on the ratio C2/C0). For L > N, larger values of 
5 < A^ do occur near the boundary at C2 = — co (32 = 0). 

Turning to L ^ A^ in region II, we have already 
pointed out that the SU(4)i state is the zero-energy GS 
for N — 3p particles at L — 3p(2p — 1) when qq and 
52 > 0, while the S0(5)i state is the zero-energy GS for 
N ^2q particles at L = 2g(g- 1) for go = 0. (For A^ = 3, 
4, these states occur at L — N). For even larger L, each 
of these model cases possesses many degenerate GS's of 
zero energy. This implies that, within our model, the 
SU(4)i and S0(5)i GS's are those found (for the param- 
eters as stated) at the critical rotation frequency w = wo, 
and so the lowest possible filling factor is 3/4 in the region 
go, 52 > 0, but is 1 when go = 0. 

At intermediate L > N values in this region, we expect 
similar physics. Thus, for C2 = co/2 > 0, {go = 0), 
the system can lower its energy by forming S0(3) singlet 
pairs of bosons. For L < 2q{q — 1) {N = 2q > 4) we do 
not have exact eigenstates, but we expect that, similar to 
Ref. Q , for V = k less than some critical value v'^> \, the 
bulk of the fluid will be in the S0(5)fe state. For C2 = 0, 
Co > 0, the preferred behavior is singlet formation via 
triples of bosons of spin zero; each such three bosons 
must be in an antisymmetric orbital state as well. For 
L < 3p{2p — 1) (A^ = 3p > 3), we do not generally have 
exact eigenstates, but we expect the bulk of the fluid 
to be the SU(4)fe states for 1/ = 3fc/4 less than another 
critical i^" > 3/4. 

We also found that the exact GS for A^ = 6, i = 6 at 
C2 = is the SU(4)2 or BTC state (|). Prompted by this, 
we have proved analytically that for N — 3p particles the 
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FIG. 2. Density profiles versus a; at y = in the LLL 
model: BTC (solid line); non-rotating BEG (broken line). 
Vertical axis is in units of the particle number N = 3p. 



SU(4)p state is an exact eigenstate of iJint with C2 — 0, 
with eigenvalue 



(9) 



We have numerically checked that for iV = 3, 6, 9, 12 
this state is the GS and believe this is true for all N — 
3p. [The energy eq. (^) is lower than the value jN{N — 
2)(27r)~'^/^co of an SU(3) multiplet of vortex states for 
spin-1 bosons, that have the same orbital wavefunction 
as the Co > vortex state in the scalar case.] 

The BTC state in eq. (||) has a clear experimental sig- 
nature, that is the density profile given in Fig. ||, whose 
form is independent of N. We caution again that one 
should probably not think of the limit = 3fc as a qH 
state, since the drop is so small. We identified the state 
as an extreme member of the SU(4)fc qH series, but it 
might be more useful to view it as a boson-triplet analog 
of the BCS paired electron states, or alternatively as an 
SU(3) analog of a skyrmion spin texture (see Ref. JlSf). 

The response of the spin-1 boson system to a rotation 
frequency u is found by minimizing Htot- In Fig- ^ we 
display the result for = 6 particles. The BTC is at 
L = 6, and the SU(4)i state is at L = 18. The degenerate 
spin multiplets listed at the steps at L < 6 each form a 
single SU(3) multiplet. The wavefunctions of these GS's 
ai L < N are uniquely determined by their SU(3) and L 
quantum numbers. 

The S0(5)i state found here will survive for some dis- 
tance off 50 — 0. For N large at fixed v, there will be sev- 
eral phases within region II, and in particular a boundary 
between the S0(5)i and SU(4)i phases that approaches 

(Jo = as w ^ Wo- 

Our constructions for C2 = directly generalize to 
an n-component rotating Bose gas with repulsive spin- 
independent (5-function interactions, implying SU(n) 
symmetry. In particular, we expect a "boson n-plet con- 
densate" at i = (n — l)A^/2, and SU(n 4- l)fc states at 
sufficiently small v — nk/{n + 1). 

To conclude, we have found several interesting states of 
matter, not containing vortices, in rotating spin-1 bosons. 

As this manuscript was prepared, two preprints ap- 
peared [|l9|,p0| that also address a rotating BEC. These 
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FIG. 3. GS values for L as a function of the rotation fre- 
quency LJ for N — 6 particles at C2 =0. (We put cjq = 1 and 
Co = (27r)^/V4-) The state at L = 6 is the BTG, while the 
L — IS state is the SU(4)i state. The spin values S are shown 
for each GS. 

papers point out the relevance of the SU(n -t- l)i states 
at large L, but do not discuss the other qH states or the 
BTC state that we consider. 
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